The cells of the sensory epithelium, and not the stria vascularis, are the main cochlear cells related to the genetic pathogenesis of age-related hearing loss.

Age-related hearing loss (ARHL) is a major health concern among the elderly population. It is hoped that increasing our understanding of its underlying pathophysiological processes will lead to the development of novel therapies. Recent genome-wide association studies (GWASs) discovered a few dozen genetic variants in association with elevated risk for ARHL. Integrated analysis of GWAS results and transcriptomics data is a powerful approach for elucidating specific cell types that are involved in disease pathogenesis. Intriguingly, recent studies that applied such bioinformatics approaches to ARHL resulted in disagreeing findings as for the key cell types that are most strongly linked to the genetic pathogenesis of ARHL. These conflicting studies pointed either to cochlear sensory epithelial or to stria vascularis cells as the cell types most prominently involved in the genetic basis of ARHL. Seeking to resolve this discrepancy, we integrated the analysis of four ARHL GWAS datasets with four independent inner-ear single-cell RNA-sequencing datasets. Our analysis clearly points to the cochlear sensory epithelial cells as the key cells for the genetic predisposition to ARHL. We also explain the limitation of the bioinformatics analysis performed by previous studies that led to missing the enrichment for ARHL GWAS signal in sensory epithelial cells. Collectively, we show that cochlear epithelial cells, not stria vascularis cells, are the main inner-ear cells related to the genetic pathogenesis of ARHL.

The predictive capacity of polygenic risk scores for disease risk is only moderately influenced by imputation panels tailored to the target population.

MOTIVATION: Polygenic risk scores (PRSs) predict individuals' genetic risk of developing complex diseases. They summarize the effect of many variants discovered in genome-wide association studies (GWASs). However, to date, large GWASs exist primarily for the European population and the quality of PRS prediction declines when applied to other ethnicities. Genetic profiling of individuals in the discovery set (on which the GWAS was performed) and target set (on which the PRS is applied) is typically done by SNP arrays that genotype a fraction of common SNPs. Therefore, a key step in GWAS analysis and PRS calculation is imputing untyped SNPs using a panel of fully sequenced individuals. The imputation results depend on the ethnic composition of the imputation panel. Imputing genotypes with a panel of individuals of the same ethnicity as the genotyped individuals typically improves imputation accuracy. However, there has been no systematic investigation into the influence of the ethnic composition of imputation panels on the accuracy of PRS predictions when applied to ethnic groups that differ from the population used in the GWAS. RESULTS: We estimated the effect of imputation of the target set on prediction accuracy of PRS when the discovery and the target sets come from different ethnic groups. We analyzed binary phenotypes on ethnically distinct sets from the UK Biobank and other resources. We generated ethnically homogenous panels, imputed the target sets, and generated PRSs. Then, we assessed the prediction accuracy obtained from each imputation panel. Our analysis indicates that using an imputation panel matched to the ethnicity of the target population yields only a marginal improvement and only under specific conditions. AVAILABILITY AND IMPLEMENTATION: The source code used for executing the analyses is this paper is available at https://github.com/Shamir-Lab/PRS-imputation-panels.

Shared and organ-specific gene-expression programs during the development of the cochlea and the superior olivary complex.

The peripheral and central auditory subsystems together form a complex sensory network that allows an organism to hear. The genetic programs of the two subsystems must therefore be tightly coordinated during development. Yet, their interactions and common expression pathways have never been systematically explored. MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression and are essential for normal development of the auditory system. We performed mRNA and small-RNA sequencing of organs from both auditory subsystems at three critical developmental timepoints (E16, P0, P16) to obtain a comprehensive and unbiased insight of their expression profiles. Our analysis reveals common and organ-specific expression patterns for differentially regulated mRNAs and miRNAs, which could be clustered with a particular selection of functions such as inner ear development, Wnt signalling, K+ transport, and axon guidance, based on gene ontology. Bioinformatics detected enrichment of predicted targets of specific miRNAs in the clusters and predicted regulatory interactions by monitoring opposite trends of expression of miRNAs and their targets. This approach identified six miRNAs as strong regulatory candidates for both subsystems. Among them was miR-96, an established critical factor for proper development in both subsystems, demonstrating the strength of our approach. We suggest that other miRNAs identified by this analysis are also common effectors of proper hearing acquirement. This first combined comprehensive analysis of the developmental program of the peripheral and central auditory systems provides important data and bioinformatics insights into the shared genetic program of the two sensory subsystems and their regulation by miRNAs.

Metformin Protects Against Noise-Induced Hearing Loss in Male Mice.

HYPOTHESIS: Metformin treatment will protect mice from noise-induced hearing loss (NIHL). BACKGROUND: We recently identified metformin as the top-ranking, Food and Drug Administration-approved drug to counter inner ear molecular changes induced by permanent threshold shift-inducing noise. This study is designed to functionally test metformin as a potential otoprotective drug against NIHL. METHODS: Male and female B6CBAF1/J mice were obtained at 7 to 8 weeks of age. A cohort of the females underwent ovariectomy to simulate menopause and eliminate the effect of ovarian-derived estrogens. At 10 weeks of age, mice underwent a permanent threshold shift-inducing noise exposure (102.5 or 105 dB SPL, 8-16 kHz, 2 h). Auditory brainstem response (ABR) thresholds were obtained at baseline, 24 h after noise exposure, and 1 week after noise exposure. Mice were administered metformin (200 mg/kg/d) or a saline control in their drinking water after the baseline ABR and for the remainder of the study. After the 1-week ABR, mice were euthanized and cochlear tissue was analyzed. RESULTS: Metformin treatment reduced the 1-week ABR threshold shift at 16 kHz ( p < 0.01; d = 1.20) and 24 kHz ( p < 0.01; d = 1.15) as well as outer hair cell loss in the 32-45.5 kHz range ( p < 0.0001; d = 2.37) in male mice. In contrast, metformin treatment did not prevent hearing loss or outer hair cell loss in the intact or ovariectomized female mice. CONCLUSIONS: Metformin exhibits sex-dependent efficacy as a therapeutic for NIHL. These data compel continued investigation into metformin's protective effects and demonstrate the importance of evaluating the therapeutic efficacy of drugs in subjects of both sexes.

SWI/SNF complexes are required for retinal pigmented epithelium differentiation and for the inhibition of cell proliferation and neural differentiation programs.

During embryonic development, tissue-specific transcription factors and chromatin remodelers function together to ensure gradual, coordinated differentiation of multiple lineages. Here, we define this regulatory interplay in the developing retinal pigmented epithelium (RPE), a neuroectodermal lineage essential for the development, function and maintenance of the adjacent retina. We present a high-resolution spatial transcriptomic atlas of the developing mouse RPE and the adjacent ocular mesenchyme obtained by geographical position sequencing (Geo-seq) of a single developmental stage of the eye that encompasses young and more mature ocular progenitors. These transcriptomic data, available online, reveal the key transcription factors and their gene regulatory networks during RPE and ocular mesenchyme differentiation. Moreover, conditional inactivation followed by Geo-seq revealed that this differentiation program is dependent on the activity of SWI/SNF complexes, shown here to control the expression and activity of RPE transcription factors and, at the same time, inhibit neural progenitor and cell proliferation genes. The findings reveal the roles of the SWI/SNF complexes in controlling the intersection between RPE and neural cell fates and the coupling of cell-cycle exit and differentiation.

Accelerated replicative senescence of ataxia-telangiectasia skin fibroblasts is retained at physiologic oxygen levels, with unique and common transcriptional patterns.

The genetic disorder, ataxia-telangiectasia (A-T), is caused by loss of the homeostatic protein kinase, ATM, and combines genome instability, tissue degeneration, cancer predisposition, and premature aging. Primary fibroblasts from A-T patients exhibit premature senescence when grown at ambient oxygen concentration (21%). Here, we show that reducing oxygen concentration to a physiological level range (3%) dramatically extends the proliferative lifespan of human A-T skin fibroblasts. However, they still undergo senescence earlier than control cells grown under the same conditions and exhibit high genome instability. Comparative RNA-seq analysis of A-T and control fibroblasts cultured at 3% oxygen followed by cluster analysis of differentially expressed genes and functional enrichment analysis, revealed distinct transcriptional dynamics in A-T fibroblasts senescing in physiological oxygen concentration. While some transcriptional patterns were similar to those observed during replicative senescence of control cells, others were unique to the senescing A-T cells. We observed in them a robust activation of interferon-stimulated genes, with undetected expression the interferon genes themselves. This finding suggests an activation of a non-canonical cGAS-STING-mediated pathway, which presumably responds to cytosolic DNA emanating from extranuclear micronuclei detected in these cells. Senescing A-T fibroblasts also exhibited a marked, intriguely complex alteration in the expression of genes associated with extracellular matrix (ECM) remodeling. Notably, many of the induced ECM genes encode senescence-associated secretory phenotype (SASP) factors known for their paracrine pro-fibrotic effects. Our data provide a molecular dimension to the segmental premature aging observed in A-T patients and its associated symptoms, which develop as the patients advance in age.

MYC Induces Immunotherapy and IFNγ Resistance Through Downregulation of JAK2.

Immunotherapy has revolutionized the treatment of advanced melanoma. Because the pathways mediating resistance to immunotherapy are largely unknown, we conducted transcriptome profiling of preimmunotherapy tumor biopsies from patients with melanoma that received PD-1 blockade or adoptive cell therapy with tumor-infiltrating lymphocytes. We identified two melanoma-intrinsic, mutually exclusive gene programs, which were controlled by IFNγ and MYC, and the association with immunotherapy outcome. MYC-overexpressing melanoma cells exhibited lower IFNγ responsiveness, which was linked with JAK2 downregulation. Luciferase activity assays, under the control of JAK2 promoter, demonstrated reduced activity in MYC-overexpressing cells, which was partly reversible upon mutagenesis of a MYC E-box binding site in the JAK2 promoter. Moreover, silencing of MYC or its cofactor MAX with siRNA increased JAK2 expression and IFNγ responsiveness of melanomas, while concomitantly enhancing the effector functions of T cells coincubated with MYC-overexpressing cells. Thus, we propose that MYC plays a pivotal role in immunotherapy resistance through downregulation of JAK2.

Transcriptional profiling of the response to starvation and fattening reveals differential regulation of autophagy genes in mammals.

Nutrient deprivation (starvation) induced by fasting and hypercaloric regimens are stress factors that can influence cell and tissue homeostasis in mammals. One of the key cellular responses to changes in nutrient availability is the cell survival pathway autophagy. While there has been much research into the protein networks regulating autophagy, less is known about the gene expression networks involved in this fundamental process. Here, we applied a network algorithm designed to analyse omics datasets, to identify sub-networks that are enriched for induced genes in response to starvation. This enabled us to identify two prominent active modules, one composed of key stress-induced transcription factors, including members of the Jun, Fos and ATF families, and the other comprising autophagosome sub-network genes, including ULK1. The results were validated in the brain, liver and muscle of fasting mice. Moreover, differential expression analysis of autophagy genes in the brain, liver and muscle of high-fat diet-exposed mice showed significant suppression of GABARAPL1 in the liver. Finally, our data provide a resource that may facilitate the future identification of regulators of autophagy.

The LHX2-OTX2 transcriptional regulatory module controls retinal pigmented epithelium differentiation and underlies genetic risk for age-related macular degeneration.

Tissue-specific transcription factors (TFs) control the transcriptome through an association with noncoding regulatory regions (cistromes). Identifying the combination of TFs that dictate specific cell fate, their specific cistromes and examining their involvement in complex human traits remain a major challenge. Here, we focus on the retinal pigmented epithelium (RPE), an essential lineage for retinal development and function and the primary tissue affected in age-related macular degeneration (AMD), a leading cause of blindness. By combining mechanistic findings in stem-cell-derived human RPE, in vivo functional studies in mice and global transcriptomic and proteomic analyses, we revealed that the key developmental TFs LHX2 and OTX2 function together in transcriptional module containing LDB1 and SWI/SNF (BAF) to regulate the RPE transcriptome. Importantly, the intersection between the identified LHX2-OTX2 cistrome with published expression quantitative trait loci, ATAC-seq data from human RPE, and AMD genome-wide association study (GWAS) data, followed by functional validation using a reporter assay, revealed a causal genetic variant that affects AMD risk by altering TRPM1 expression in the RPE through modulation of LHX2 transcriptional activity on its promoter. Taken together, the reported cistrome of LHX2 and OTX2, the identified downstream genes and interacting co-factors reveal the RPE transcription module and uncover a causal regulatory risk single-nucleotide polymorphism (SNP) in the multifactorial common blinding disease AMD.

Alternative cleavage and polyadenylation generates downstream uncapped RNA isoforms with translation potential.

The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Incorporating regulatory interactions into gene-set analyses for GWAS data: A controlled analysis with the MAGMA tool.

To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such "augmented analyses", and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries.

Gene architecture directs splicing outcome in separate nuclear spatial regions.

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

The DOMINO web-server for active module identification analysis.

MOTIVATION: Active module identification (AMI) is an essential step in many omics analyses. Such algorithms receive a gene network and a gene activity profile as input and report subnetworks that show significant over-representation of accrued activity signal ('active modules'). Such modules can point out key molecular processes in the analyzed biological conditions. RESULTS: We recently introduced a novel AMI algorithm called DOMINO and demonstrated that it detects active modules that capture biological signals with markedly improved rate of empirical validation. Here, we provide an online server that executes DOMINO, making it more accessible and user-friendly. To help the interpretation of solutions, the server provides GO enrichment analysis, module visualizations and accessible output formats for customized downstream analysis. It also enables running DOMINO with various gene identifiers of different organisms. AVAILABILITY AND IMPLEMENTATION: The server is available at http://domino.cs.tau.ac.il. Its codebase is available at https://github.com/Shamir-Lab.

CT-FOCS: a novel method for inferring cell type-specific enhancer-promoter maps.

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer-promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome.

Lineage-tracing and translatomic analysis of damage-inducible mitotic cochlear progenitors identifies candidate genes regulating regeneration.

Cochlear supporting cells (SCs) are glia-like cells critical for hearing function. In the neonatal cochlea, the greater epithelial ridge (GER) is a mitotically quiescent and transient organ, which has been shown to nonmitotically regenerate SCs. Here, we ablated Lgr5+ SCs using Lgr5-DTR mice and found mitotic regeneration of SCs by GER cells in vivo. With lineage tracing, we show that the GER houses progenitor cells that robustly divide and migrate into the organ of Corti to replenish ablated SCs. Regenerated SCs display coordinated calcium transients, markers of the SC subtype inner phalangeal cells, and survive in the mature cochlea. Via RiboTag, RNA-sequencing, and gene clustering algorithms, we reveal 11 distinct gene clusters comprising markers of the quiescent and damaged GER, and damage-responsive genes driving cell migration and mitotic regeneration. Together, our study characterizes GER cells as mitotic progenitors with regenerative potential and unveils their quiescent and damaged translatomes.

A cell-type-specific atlas of the inner ear transcriptional response to acoustic trauma.

Noise-induced hearing loss (NIHL) results from a complex interplay of damage to the sensory cells of the inner ear, dysfunction of its lateral wall, axonal retraction of type 1C spiral ganglion neurons, and activation of the immune response. We use RiboTag and single-cell RNA sequencing to survey the cell-type-specific molecular landscape of the mouse inner ear before and after noise trauma. We identify induction of the transcription factors STAT3 and IRF7 and immune-related genes across all cell-types. Yet, cell-type-specific transcriptomic changes dominate the response. The ATF3/ATF4 stress-response pathway is robustly induced in the type 1A noise-resilient neurons, potassium transport genes are downregulated in the lateral wall, mRNA metabolism genes are downregulated in outer hair cells, and deafness-associated genes are downregulated in most cell types. This transcriptomic resource is available via the Gene Expression Analysis Resource (gEAR; https://umgear.org/NIHL) and provides a blueprint for the rational development of drugs to prevent and treat NIHL.

A CRISPR knockout screen reveals new regulators of canonical Wnt signaling.

The Wnt signaling pathways play fundamental roles during both development and adult homeostasis. Aberrant activation of the canonical Wnt signal transduction pathway is involved in many diseases including cancer, and is especially implicated in the development and progression of colorectal cancer. Although extensively studied, new genes, mechanisms and regulatory modulators involved in Wnt signaling activation or silencing are still being discovered. Here we applied a genome-scale CRISPR-Cas9 knockout (KO) screen based on Wnt signaling induced cell survival to reveal new inhibitors of the oncogenic, canonical Wnt pathway. We have identified several potential Wnt signaling inhibitors and have characterized the effects of the initiation factor DExH-box protein 29 (DHX29) on the Wnt cascade. We show that KO of DHX29 activates the Wnt pathway leading to upregulation of the Wnt target gene cyclin-D1, while overexpression of DHX29 inhibits the pathway. Together, our data indicate that DHX29 may function as a new canonical Wnt signaling tumor suppressor and demonstrates that this screening approach can be used as a strategy for rapid identification of novel Wnt signaling modulators.

Genetic mapping of developmental trajectories for complex traits and diseases.

Genome-wide association studies (GWAS) have identified numerous common genetic variants associated with complex human traits and diseases. However, the translation of GWAS discoveries into biological and clinical insights is highly challenging. In this study, we present a novel bioinformatics approach for enhancing the functional interpretation of GWAS signals, based on their integration with single-cell (sc)RNA-seq datasets that examine developmental processes. Our approach performs three tasks: (1) Identification of links between cell differentiation trajectories and traits; (2) Elucidation of biological processes and molecular pathways that underlie such trajectory-trait links; and (3) Prioritization of target genes that carry the links between trajectories, pathways and traits. We applied our method to a set of 11 traits of various pathologies, and 12 scRNA-seq datasets of diverse developmental processes, and it readily detected well-established biological connections, including those between the maturation of cortical inhibitory interneurons and schizophrenia, hepatocytes and cholesterol levels, and pancreatic beta-islet cells and type-2 diabetes. For each of these associations, our method pinpointed top candidate genes that are strongly associated with both the kinetics of the differentiation trajectory and the disease's genetic risk. By the identification of trajectory-disease links, molecular pathways that underlie them and prioritizing candidate risk genes, our method improves the understanding of the etiology of complex diseases, and thus holds promise for enhancing rational drug development that is aimed at targeting specific biological processes that mediate the genetic predisposition to diseases.